Double-Occupancy Errors, Adiabaticity, and Entanglement 
of Spin-Qubits in Quantum Dots 
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Quantum gates that temporarily increase singlet-triplet 
splitting in order to swap electronic spins in coupled quantum 
dots, lead inevitably to a finite double-occupancy probability 
for both dots. By solving the time-dependent Schrodinger 
equation for a coupled dot model, we demonstrate that 
this does not necessarily lead to quantum computation er- 
rors. Instead, the coupled dot ground state evolves quasi- 
adiabatically for typical system parameters so that the 
double-occupancy probability at the completion of swapping 
is negligibly small. We introduce a measure of entanglement 
which explicitly takes into account the possibilty of double 
occupancies and provides a necessary and sufficient criterion 
for entangled states. 

PACS numbers: 85.30.Vw, 85.30.Wx, 03.67.-Lx, 03.67.-a, 
73.23.Hk 



I. INTRODUCTION 

In the past several years there has been a great deal 
of interest in possible physical realizations of quantum 
computing bits and operations B. Among the various 
proposals, solid state systems are particularly attractive 
since they are more easily integrated into large quantum 
networks. In particular, semiconductor nanostructures 
which use the spin degree of freedom of the electrons 0] 
(rather than their charge) for information processing are 
of special interest since they can take advantage of the 
comparatively long spin coherence times in such materi- 
als |-|. 

A key challenge is the construction of systems com- 
posed of two coupled quantum dots which can be coupled 
to perform swap operations Usw, he. unitary two-qubit 
operations which interchange the spin states (qubits) of 
the electrons on the two dots |p|,p|-pOl . By combining the 

1 /2 

"square root" Ug-^ of such a swap with other isolated- 
qubit manipulations one can construct a quantum XOR 
gate. A quantum XOR gate, along with isolated-qubit 
operations, has been shown to be sufficient for the im- 
plementation of any quantum algorithm ||ll|. Hence a 
practical and reliable realization of a swap gate would be 
an important step towards the fabrication of a solid state 
quantum computer. 

The swap operation of electron spin states in a double 
dot system can be realized in principle by turning on a 
time-dependent exchange coupling between the spins as a 



"source of entanglement" . In practice the exchange inter- 
action is provided by singlet-triplet splitting in a double- 
dot, which is always accompanied by a finite inter-dot 
electron tunneling amplitude. 

In a recent work Das Sarma and coworkers |§,|l^ em- 
phasized that exchange interactions in the range of inter- 
est are accompanied by a substantial probability, during 
the swap operation, that both electrons will be on the 
same dot. In this paper we demonstrate that, contrary 
to naive expectations, these virtual double-occupancies 
will not, under circumstances typically envisioned, lead 
to an important increase in quantum computing errors. 
Double occupancy is not a fatal problem for quantum dot 
based quantum computing with spins. The occurrence of 
double occupancies during the swapping process does not 
lead to processing errors, provided that the double occu- 
pancies are sufficiently suppressed when the swapping of 
spin states is completed. The principle purpose of the 
present paper is to illustrate this basic feature within the 
Hund-MuUiken description of a quantum dot hydrogen 
molecule. We will see that, in a system of identical dots, 
the time evolution of this system can be reduced to the 
problem of a pseudospin-i in a time-dependent pseudo- 
magnetic field. In particular, the question of whether 
double occupancies are a severe obstacle for swap oper- 
ations in the quantum dot system is equivalent to the 
question of how close the pseudospin dynamics is to its 
adiabatic limit. Simple numerical studies presented in 
Sec. IV show that the pseudospin has an approximately 
adiabatic time evolution for a ramarkably broad range 
of coupling ramp times. It turns out that this behav- 
ior holds even if the inversion symmetry along the x-axis 
connecting the dots is broken (e.g. in the presence of an 
electric field). 

A secondary purpose of this work is to introduce a 
coordinate-independent measure of entanglement appro- 
priate for the Hilbcrt space of the above system. This 
quantity provides a necessary and sufficient criterion for 
the entanglement of quantum states. It differs from other 
entanglement criteria proposed in the literature |13| in so 
far as it takes into account states with double occupan- 
cies. This generalizes the typical situation of Einstein- 
Podolsky-Rosen experiments. We expect this measure of 
entanglement to be useful in the theoretical study of cou- 
pled quantum dots (or similar quantum-confined nanos- 
tructures), independent of the particular model consid- 
ered here. 



II. THE MODEL 

We consider a system of two electrons in two laterally 
coupled quantum dots. The experimental motivation for 
the model described below has been discussed elsewhere 
|^-||. Here we just summarize its basic features. 

The Hamiltonian is given hy H — T + C, where C de- 
notes the Coulomb repulsion between the electrons, and 



T = E^ 



, hi is the one-particle part with 



'^i = T^(p^ + -An) ) + V{n) . 

Zm \ c ' 



(1) 



The one-particle Hamiltonian hi describes electron dy- 
namics confined to the xy-plane in a perpendicular mag- 
netc field B. The effective mass m is a material-dependent 
parameter. The coupling of the dots (which includes tun- 
neling) is modeled by a quartic potential 



V{x,y)^ 



1 

4^ 



'x^-a^f 



V 



(2) 



which separates into two harmonic wells of frequency ojq 
(one for each dot) in the limit 2a 3> 2ao, where a is half 
the distance between the dots and a^ = yJh/mujQ is the 
effective Bohr radius of a dot. 

Following Burkard et al. M we employ the Hund- 
MuUiken method of molecular orbits to describe the low- 
lying spectrum of our system. This approach concen- 
trates on the lowest orbital states in each dot and is an 
extension of the Heitler-London method also discussed 
in H]. The Hund-MuUiken approach accounts for dou- 
ble occupancies and is therefore suited to investigate the 
questions at issue here. 

In the usual symmetric gauge A = B{~y, x,0)/2 the 
Fock-Darwin ground state of a single dot with harmonic 
confinement centered around r — (±a,0, 0) reads 



muj fmuj f 2 



I a 

•expi±2yp- 



(3) 



where Is — \/hc/eB is the magnetic length, and the 
frequency uj is given by to^ = cjq + u}\ where 



^L 



eB /2mc is the usual Larmor frequency. From these non- 
orthogonal one-particle states we construct the orthonor- 
malized states \A) and \B) with wavefunctions 



[f\A) = 



1 



[ip+a - g^-a) , 



(4) 



v/1 - 2Sg - g2 

{■AB) = ^ ^^ ((p-a~9f+a), (5) 



v/1 - 2Sg - .g2 
with S being the overlap between the states 



and 



Vl — S-^)/S. For appropriate values of system 
parameters such as the interdot distance and the exter- 
nal magnetic field, the overlap S becomes exponentially 



small [^. In this limit an electron in one of the states \A), 
\B) is predominantly localized around r = (±a, 0, 0). In 
the following we consider this case and use these states 
as basis states to define qubits, i.e. qubits are realized 
by the spin state of an electron in cither orbital \A), or 
orbital \B). 

An appropriate basis set for the six-dimensional two- 
particle Hilbert space is given (using standard notation) 
by the three spin singlets 



l^i> = 



1 



V2 

and the triplet multiplet, 
1 



^AI^Al 



^AI^Al 



^B-[^Bi 






|0), 
|0), 
|0), 



1^/^/2 ['^A]'^Bi + ^Ai^B] 



|0), 



(6) 

(7) 
(8) 

(9) 
(10) 

(11) 



The three triplet states are degenerate (typically we can 
ignore possible Zeeman splittings p]) and have the com- 
mon eigenvalue. 



er = 2£ + VL , 
where we have defined 



e = {A\h\A) = {B\h\B) 



and 



V- = (T"|c|r") 



V+ = {Si\C\Si) 



(12) 



(13) 



(14) 



An important further observation is that, as a conse- 
quence of inversion symmetry along the axis connecting 
the dots, the Hamiltonian does not have any non-zero 
matrix elements between the singlet state jS'a) and other 
states. Hence, jS'a) is, independently of the system pa- 
rameters, an eigenstate. The eigenvalues of the triplet 
and jS'a), however, do depend on system parameters. The 
Hamiltonian acting on the remaining space spanned by 
IS*!) and 1 52) can be written as 



n = 2e+-UH + V+- 



( UhI2 2tH 
{ 2tH -Uh/2 



where 



and 



tH = {A\h\B) = {B\h\A) 



Uh = {S2\C\S2)^V+. 



(15) 



(16) 



(17) 



The nontrivial part of ( |15D is a simple Hubbard Hamil- 
tonian and can be identified as the Hamiltonian of a 
pseudospin-i-object in a pseudomagnetic field having 
a component Ur in the z-direction and Atn in the x- 
direction of pseudospin space. [Note that this pseudospin 
is not related to the spin degree of freedom which pro- 
vides the qubit!] The space spanned by \Si) and IS'2) 
contains the ground state of the system. The basis states 
themselves are eigenstates only in the case of vanishing 
tunneling amplitude tn where \Si) is the ground state. 
In all other cases, the ground state has an admixture of 
double occupied states contained in IS'2)- The energy gap 
between the triplet and the singlet ground state is 

eT-es^ = V+-V^-^ + ^Uf, + Wtj, . (18) 

A swap operation in the present system is a unitary 
transformation which turns a state having the qubits in 
different states, say. 



1 



c%ct,\0) = ^{\T') + \S,)) 



(19) 



into a state where the contents of the qubits is inter- 
changed, 



1 



= ^ {\T') \Si)) 



(20) 



These two states are eigenstates in the case V+ = V- and 
tn = for which the singlet-triplet splitting vanishes. 

As discussed in references ||,^, swapping may be 
achieved by the action of a gate that lowers the poten- 
tial barrier between the quantum dots. This leads to 
exponentially larger values for both V+ — V- and tn- It 
is adequate for our purposes to consider a model where 
V+ = V- (consistent with the above limit of small overlap 
iS"), and the singlet-triplet splitting results entirely from 
tn- If the duration and amplitude of a tunneling pulse 
is adjusted appropriately, the relative phase between the 
singlet and the triplet state involved picks up a shift of 
TT, and a swapping operation is performed. 

As pointed out in |10| a finite tunneling amplitude nec- 
essarily leads to a finite probability for double occupan- 
cies of qubit states. If double occupancy errors occurs to 
any sizable extent as a result of the swapping process, any 
quantum computation based on this hardware is likely to 
fail. However, if the double occupancies are sufficiently 
rare after the swapping process, errors in the quantum 
computation can likely be corrected dynamically. An im- 
portant observation is that the double-occupancy proba- 
bility after the swap vanishes in the adiabatic limit, i.e. if 
the ramp time r of the quantum gate is such that h/r is 
much larger than the pseudospin splitting y^Un + 16i|^. 
This follows since the non-adiabatic effects can arise only 
from the states \Si) and IS'2), which have a non-trivial 
time evolution when the tunneling amplitude tn is time- 
dependent. Thus, the question of whether double oc- 
cupancies are problematic for swap operations in the 



present system is reduced to the question of how close the 
motion of a spin-i-object in a time-dependent magnetic 
field is to its adiabatic limit. This will be investigated 
further in Sec. |^. 

The reduction of the dynamics to the time evolution 
of a two-level-system relies on the fact that the system 
has inversion symmetry along the x-axis in real space 
connecting the dots. This symmetry can be broken if 
odd powers of the particle coordinates Xi are added to 
the Hamiltonian (|l|) like for example the potential of a 
homogeneous electric field. However, the only addional 
matrix element due to such terms in the Hamiltonian 
occurs in the subspace of double occupied states between 
the singlets IS2) and IS'3). Thus, in the presence of an 
electric field £ — —eE'Y^^Xi the Hamitonian acting on 
the singlet subspace spanned by |Si), IS2), IS3) reads 



71 = 26 +-UH + V+ 



( Uh/2 2tH 

2tH -Uh/2 F 

\ F -Uh/2 + 2X 



(21) 



with the real matrix element F = (S2|f IS3) and 

2X = {S2\C\S2) - (S3ICIS3) = 2{A\{A\C\B)\B) . (22) 

With a finite matrix element F the dynamics of the sys- 
tem is slightly more complicated, but also in this case 
the only coupling of the two-qubit states ([l^) and ( [20| ) 
to the subspace of double occupied states is provided by 
the tunneling amplitude tn- Therefore, with respect to 
the adiabaticity of the swapping process, the situation 
can be expected to be not very different from the one 
with inversion symmetry between the dots. This will be 



verified in Sec. IV 



So far we have not considered a possible Zeeman cou- 
pling to the electron spin. This would not change the 
situation essentially since all states involved in the swap- 
ping process (|r"), |Si), IS'2), and eventually IS'3)) have 
the total spin quantum number S^ = 0. 



III. ENTANGLED STATES 

Before analysing further the possibility of performing 
swap operations in the above system, let us introduce an 
appropriate measure for the entanglement of its quantum 
states. Consider a system of two fermions living in a four- 
dimensional one-particle space. A general state vector in 
this six-dimensional Hilbert space can be written as 



\w) = WafcC+C+|0) 



(23) 



where a summation convention is understood for re- 
peated latin indices a, 6, ... G {1, 2, 3, 4} running over the 
orthonormalised one-particle states. The coefficient ma- 
trix w can be assumed to be antisymmetric, Wab = —Wba- 
The normalisation condition reads 



{w\w) =1 <^ Tt{ww) = — - . 



(24) 



where the bar denotes complex conjugation. A two- 
particle state of the form (|2|) is in general entangled, 
i.e. cannot be written as a single Slater determinant. 
\w) is non-entangled, i.e. a single Slater determinant, if 
w has the form 



Wab = 77 [ZaZb 



%^a) 



(25) 



for two orthonormal spinors z^ and z^, z^^zl = 5^K We 
note that for a given non-entangled state \'w) the choice 
of spinors zi and Z2 is by no means unique since any 
SU(2) transformation among these two occupied one- 
particle states leads to the same two-particle state vector 
\w). Hence, for a given non-entangled state \'w) there is a 
three-dimensional manifold of spinors fulfilling equation 

We define the dual matrix w of w by 



1 



Wab 



abed - 

£ Wed 



(26) 



with e"'^'^'^ denoting the totally antisymmetric unit tensor 
in four dimensions. The scalar product of a state \w) with 
its dual state \'w) can be written as 

{w\w) = £'''""^WabWcd = 8 (^12^34 + W13W42 + W14W23) ■ 

(27) 



This cyclic sum vanishs identically if w has the form (25) 
Hence the vanishing of 



r]{w) :— \{w\w) 



(28) 



is a necessary condition for \w) being a single Slater de- 
terminant. Moreover, in the Appendix it is shown that 
r]{w) = is actually also a suffcient condition for |w) be- 
ing non-entangled. Thus, rj defines a measure of entan- 
glement which is exactly zero for non-entangled states. 
r]{w) 7^ is therefore a necessary and sufficient condition 
for entanglement of quantum states. Maximally entan- 
gled states are characterized by the fact that they are 
collinear with their dual states, rj^w) — I. As simple 
examples, we consider the basis states used in the pre- 
ceding section: the states \T~^) and \T^) are single Slater 
determinants, while all other basis states are maximally 
entangled. 

The matrix w transforms under a unitary transforma- 
tion of the one-particle space. 



as 



■ UC+U+ = Ubact 



w 1-^ UwU 



(29) 



(30) 



where U"^ is the transpose (not the adjoint) of U. It 
is straightforward to see that rj is invariant under such 



transformations, and the determinant of w remains the 
same up to a possible phase factor. Thus, the entangle- 
ment of a state {w) quantified by 77 does not depend on 
the basis chosen for the one-particle space, which is of 
course a necessary requirement for a measure of entan- 
glement. 

The dualisation of a state can be identified as a 
particle- hole-transformation, 

Up-hctU+_^ ^ Ca , l(p-h\0)^C+C+C+C+\0), (31) 

followed by complex conjugation. In fact. 



Up-h\w) = -\w) 



(32) 



We note that the complex conjugations in (Eq) and (p2) 
are unimportant for states \w) such that 77(1/7) — 0, since 
a single Slater determinant is always orthogonal to its 
particle-hole conjugate, irrespective of a possible phase 
transformation of a prefactor. However, the complex 
conjugation in the definition ( P6|) is essential to the suf- 
ficiency of the above condition. As an example consider 
a state \w) with W12 = W34 = 1/4, W13 = u'24 = i/'^ 
and wi4^ — W2Z = 0. This is clearly a maximally entan- 
gled state, rj{w) — 1, while its scalar product with the 
complex conjugate of its dual state is {'w\w) = 0. 

We also mention the following identity for the deter- 
minant of w. 



det w = 



(33) 



Hence it follows that also | det w\ could be used as a mea- 
sure of entaglement. Equation ( p3[ ) is important for the 
proof of the sufficiency of our criterion for non-entangled 
states, as explained in detail in the Appendix. 

A convenient choice to make contact between the gen- 
eral state labels a, 6,... G {1,2,3,4} used here and 
the basis states of the preceding section is given by 
(1,2,3,4) = {A ],A i,B ],B [). With this conven- 
tion, a state vector spanned by IS'2) and IS'3) only has 
W12 and i(;34 as its only independent non-zero coefficients 
in w. Such a state lies fully in the subspace of double 
occupancies, and its entanglement is purely due to the 
orbital degrees of freedom: 



r]orb = 8|'Wl2W34| 



(34) 



On the other hand, a state spanned by \Si) and jT") 
has no double occupancies and is entangled purely with 
respect to the spin degrees of freedom. 



Vspin = ^\wiaW2z\ 



(35) 



For a general state vector, both kinds of entanglement 
(orbital and spin) contribute to ri{'w). 







IV. RESULTS FOR THE SWAPPING PROCESS 

We now continue with our investigation of the dynam- 
ics of the double quantum dot qubit swapping process 
generated by a time-dependent tunnehng amphtude. 

Let us first consider the case of inversion symmetry 
along the axis connecting the dots. As explained in Sec. 

this problem can be reduced essentially to the time 
evolution of a pseudospin-i -object in a magnetic field 
having a time-dependent component in the x-direction 
of the pseudospin space. In the course of swapping, the 
triplet contribution to the incoming state (M) will just 
pick up a phase factor according to its constant eigen- 
value, while the singlet contribution will mix with the 
other singlet 152). Therefore, a finite probability for dou- 
ble occupancies will necessarily arise during the swapping 
process. However, if these amplitudes can be suppressed 
sufficiently when the swapping is complete (as in the adi- 
abatic limit) , errors in the quantum computation can be 
avoided. Thus we are left with the question of how close 
the dynamics of our formal spin- ^ -object is to its adi- 
abatic limit. We note that, in the adiabatic limit, no 
Berry phase occurs in the time evolution of the singlet 
states, since the motion of the formal spin is restricted to 
a plane. Hence the solid angle encircled in a round trip 
is strictly zero. 

The integration of the Schrodinger equation for our 
time-dependent two-level problem is in general non- 
elementary. However, there is a considerable body of lit- 
erature, starting with early work by Landau M], Zener 
||14| , and Rosen and Zener |1^, where particular cases 
of this problem were reduced to well-known differential 
equation of mathematical physics such as the hyperge- 
ometric equation. This work was reviewed and gener- 
alized very recently in |16| . However, such an approach 
still works only for special time-dependent Hamiltonians, 
i.e., in the present context, only for special shapes of the 
tunneling pulse iff(i), and many quantities of interest are 
given by complicated non-elementary expressions which 
require numerical evaluation. For this reasons, and for 
the sake of brevity of our presentation, we shall resort 
to numerical integrations of the Schrodinger equation. 
From such studies we will see that the range of adia- 
baticity is remarkably large. Our numerical findings will 
be corroborated and made physically plausible by well- 
known analytical results for Landau-Zener-type transi- 
tions in simplified cases. 

To be specific, we consider a time-dependent tunneling 
of the form 



tnit) 



1 



sh(t/r 



(36) 



cosh(T/(2r)) 



This is a tunneling pulse which is switched on and off 
exponentially with a characteristic time t. It has a du- 
ration of T and an amplitude given by A (for T ^ t). 
Therefore this form is flexible enough to describe the es- 
sential features of a pulse. The exponential switching 



is motivated by the exponential-like dependence of the 
tunneling matrix element on external parameters g] . 

A typical situation is shown in figure |l| for a switch- 
ing time of T = Ah/Un, an amplitude of A = Uh/8 and 
the duration T adjusted to enable single swap operation. 
The figure shows the results of a numerical integration 
of the time-dependent Schrodinger equation using the 
fourth order Runge-Kutta scheme. The time-dependent 
tunneling amplitude tnit) is plotted (in units Uh) as a 
dotted line. The square amplitude of the incoming state 
(O) and the outgoing state (120) are shown as thick lines. 
The square amplitudes of the singlets l^i) and 15*2) are 
denoted by |(/3ip and |(/J2p, respectively, and plotted as 
long-dashed hues. The probability of double occupan- 
cies is given by |iy92p- As one can see from the figure, 
this quanity is finite during the swapping process but 
strongly suppressed afterwards. The measure of entan- 
glement ri{t) is also shown in the figure. It is zero for the 
non-entangled incoming and outgoing state, and achieves 
its maximum value of almost unity in the middle of the 
process. This quantifies and shows explicitly the entan- 
glement of the quantum state during the swapping pro- 
cess. 

The probability |(^2p for double occupancy after 
switching off the tunneling depends on the switching time 
T, the amplitude A and also on the duration T of the tun- 
neling pulse, i.e. on the exact time when the switching off 
sets in. However, our numerics suggest that there is an 
upper bound for |(^2p at given r and A. In the above ex- 
ample the double occupancy probability after the swap- 
ping process is smaller than 10^^", which is a very tiny 
value. A typical order of magnitude for the double oc- 
cupancy probability is 10^® for amplitudes A < Uh and 
switching times r > ATi/Uh- In fact, also larger values of 
A (being still comparable with Uh) can be possible, lead- 
ing to double occupancy probabilities of the same order, 
while this probability significantly increases if r becomes 
smaller than Ah/Un- Thus, this value characterizes the 
region where the motion of the system is close to its adi- 
abatic limit and is remarkably small on the natural time 
scale of the system given by h/Un, while adiabatic be- 
havior is in general expected for a particularly slow time 
evolution. 

This large range of quasi-adiabatic behavior can be 
understood qualitatively by considering a simplified situ- 
ation where the tunneling is switched on and off linearly 
in time and is constant otherwise. Then, non-adiabtic ef- 
fects can occur only during the sharply defined switching 
processes. For simplicity, we consider the first switch- 
ing process only where the tunneling has the time de- 
pendence tn = (A/t)^, t g [0,t]. To enable analyt- 
ical progress let us further assume t G [—00,00], which 
should lead to an upper bound for the probability of non- 
adiabatic transitions due to the switching. This problem 
was considered a long time ago by Landau p3[ and by 
Zener [Q . The result of reference ||lj] for the probability 
of non-adiabatic transitions reads 



Pr 



nad 



with an adiabaticity parameter 



r/2 



8 n(A/r) 



(37) 



(38) 



We see that the probability for non-adiabatic transi- 
tions is exponentially suppressed with increasing switch- 
ing time T. This exponential dependence explains quali- 
tatively the above observation of a large range of quasi- 
adiabatic bahavior. To obtain an estimate for a nonlinear 
switching one may replace the ratio (A/r) in the denom- 
inator of ( |3q ) by the maximum time derivative of the 
tunnehng t//(i) (giving a = TTUfj/3h{A/T) for the pulse 

A similar exponential dependence of the probability for 
non-adiabatic transitions on the switching time r was 
also found analytically by Rosen and Zener |15[ for a 
particular two-parametric pulse of the form 



tnit) = A/cosh(t/T). 



(39) 



In this case non-adiabatic transitions occur with a prob- 
ability 



Pr 



nad 



sin^ (At/(2^)) / cosh^ {UHT/{2h)) 



(40) 



To illustrate the behavior in the strongly non-adiabatic 
case we have plotted in figure \'^i\'^ and |(^2p for the 
same situation as in figure |l| but with a four times smaller 
ramp time of only r = Tl/Uh- In this case small oscilla- 
tions occur in the time-evolution of these two quantities 
during the tunneling pulse, which can be understood in 
terms of the eigenspectrum at a given tunneling tn — A- 
Additionally, a sizable double occupancy probability of 
about 0.005 remains after the pulse, as shown in the in- 
set. 

Figure shows a square root of a swap, which is ob- 
tained from the situation of figure |l| by halfing the du- 
ration T of the tunneling pulse. The resulting state is a 
fully entangled complex linear combination of the states 
\Si) and |T''), or, equivalently, of the incoming state ( pj| ) 
and the outgoing state (20) of the full swap. Again, the 
weight of the doubly occupied state 15*2) is strongly sup- 
pressed after the tunneling pulse. As a consequence, Eq. 
(p3) implies that rjorb = after completion of switching, 
while rj = rjspin = 8|wi4'u;23| = 1. This shows that the 
entanglement of the two electrons is entirely in the spin 
(and not in the orbital) degrees of freedom after switch- 
ing. 

Let us finally consider swapping processes when the 
inversion symmetry along the axis connecting the dots is 
broken. Such processes are governed by the Hamiltonian 
( pi| ) in the presence of a finite matrix element F . Our 
numerical results are in this case qualitatively the same 
as before with the admissible switching times t slightly 
growing with increasing F . In figure B we illustrate our 
findings for a comparatively large off-diagonal element 



F — QAUh- The additional Coulomb matrix element is 
X — 0.2Uh, and the parameters of the tunneling pulse 
are r = SH/Uh and A = Uh/8 with a duration T appro- 
priate for a single swapping. As a result, a clean swap- 
ping operation can be performed also in the absence of 
inversion symmetry. 

We note that the Hund-MuUiken scheme used here is 
restricted to the low-energy sector where only the lowest 
single-particle energy levels (with typical spacings 5e) are 
kept. For this scheme to be valid also in a switching 
process, we need to require that time-dependent changes 
must be performed adiabatically also with respect to the 
time scale set by h/Se, i.e. we need r > H/Se |^. On 
the other hand, to suppress double occupancy errors we 
have seen that the adiabaticity parameter a of Eq. (|3q ) 
must exceed one, implying that r > 8hA/{TTU]f). Thus, 
the adiabaticity condition for switching becomes more 
generally. 



T > T„ 



sn A 



(41) 



There are now two particular cases we can distinguish. 
First, if the effective Coulomb charging energy exceeds 
the level spacing, i.e. Uh > Se, we obtain Tmin — h/Se, 
since for consistency we have A < Se. Thus, when the 
switching is adiabatic with respect to the scale set by 
Se, errors due to double occupancy are automatically ex- 
cluded. In the second case with Uh < VA(5e < Se ("ul- 
trasmall quantum dots"), we obtain Tmin = 8?iA/(7rt/^), 
which means that the overall condition for adiabaticity is 



determined by the no-double occupancy criterion. Usin, 

)r GaAs quantum dots [| 
is of the order of 50 ps 



typical material parameters for GaAs quantum dots [0 , 
we can estimate [pf that t„ 



V. CONCLUSIONS 

We have studied a double quantum dot system as a 
quantum gate swapping the electronic spin states on the 
two dots. Within Hund-MuUiken approach, the dynam- 
ics of such a system having inversion symmetry along 
the axis connecting the dots reduces to the problem of 
a pseudospin-2-object in a time-dependent pseudomag- 
netic field. By solving the time-dependent Schrodinger 
equation numerically we demonstrate the possibility of 
performing swap operations and investigate the role of 
double ocupancies of the dots. These double occupancies 
are found to be (exponentially) strongly reduced, as a 
result of the swapping process, for a large range of sys- 
tem parameters and are therefore not a principle obstacle 
for quantum computation in such systems. Further nu- 
merical studies show that this situation is not altered 
qualitatively when the inversion symmetry is broken. 

Moreover, we have introduced an appropriate measure 
of entanglement which takes explicitly into account the 
possibility of double occupancies. This quantity allows to 
quantify the entanglement of the quantum state during 



a gate operation and provides a necessary and sufficient 
condition for entangled states. Hence we expect this mea- 
sure of entanglement to be useful in general in the study 
of quantum information phenomena in systems such as 
(real or artificial) diatomic molecules, or other quantum- 
confined two-site structures. 
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Since both terms in the imaginary part are equal by def- 
inition and det w = implies {v\v) = 0, it holds 



(u\u) = 



dot M = 



and 



{u\v) — {v\u) — . 



(A4) 



(A5) 



From (i) one concludes that both \u) and \v) are single 
Slater determinants. Thus there are spinors x^, x^ and 
2/^, y^ with 



Uah — 2 V^a'^h ^ ^h^a) 



1/12 1 2\ 

Vab = 2 [yaVb ~ VbVa) 



Moreover, equation (A5) implies that 



pQ&cd 1 2 1 2 _ r) 



(A6) 



(A7) 



APPENDIX A: 

Here we give the proof that ri{w) = is indeed a suff- 
cient condition for \w) being a single Slater determinant 
state. The proof consists of two steps. 

(i) Let vf be purely real. Since rj(w) = implies 
deiw = (cf. eq. (p3|)), w has at least one zero eigen- 
value. Because w is real and antisymmetric its eigenval- 
ues are purely imaginary (if not zero) and occur in pairs 
of complex conjugates. Therefore, at least two of the four 
eigenvalues of w are zero. It follows from standard argu- 
ments (similar to those for real and symmetric matrices) 
that these two zero eigenvalues correspond to two real 
eigenvectors being orthogonal onto each other. It follows 
that there is a real and orthogonal one-particle transfor- 
mation U so that, say, the first two rows and columns 
of the resulting matrix UwU'^ are zero. Hence, the one- 
particle states with labels a = 1,2 (in this new basis) 
are strictly empty, and the two electrons occupy the re- 
maining two states. Thus, {w) is clearly a single Slater 
determinant. 

(ii) General case: w complex. By a one-particle 
transformation U with 



U 



diag(e*'^\e*'^%e''^%e''^'') 



(Al) 



one can adjust the phases in w' = UwU'^ in a manner 
that, say, w[2, w'^^, w'^ are real. Denoting the real and 
imaginary part of w' by 

w' ^ u + iv (A2) 

it follows that det v = 0. Consider now the (unnor- 
malized) states \u) and |f). If one of these states van- 
ishes the assertion is already proved in (i), thus assume 
|u) 7^ 7^ \v). The condition rj{w) = reads 



Thus, the 4x4-matrix having these four spinors as its 
rows or columns has a vanishing determinant. Therefore 
these spinors are linearly dependent. Without loss of 
generality, consider the case 



Py' 



ly 



(A8) 



where the complex coefficients (3 and 7 are not both zero 
since otherwise u = Q. Let, again without loss of gener- 
ality, (3 be nonzero. Then the spinors 



z'^ = (3y^ + 72/2 
solve the problem, i.e. 



x^ + -/ 



Uab + iVab = 2 i^-l^-b - 4^'l) 



(A9) 



(AlO) 



z^ and 2^ are both nonzero and not coUincar to each other 
since otherwise w' — 0. Thus, up to an unimportant 
orthonormalisation, these two spinors define one-particle 
states which allow to express \w') (and consequently \w)) 
as a single Slater determinant. 



\u) - {v\v) + i {{u\v) + {v\u)) = 



(A3) 
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FIG. 1. A swap process as a function of time. The tunnel- 
ing amphtude tH{t) is plotted (in units of Uh) as a dotted 
line. The square amplitude of the incoming state ( |l9[ ) and 
the outgoing state (EOl) are shown as thick hues. The square 
amplitudes of the singlets jSi) and j^) are denoted by \ipi\ 
and \^2\^ , respectively, and plotted as long-dashed lines. The 
measure of entanglement ri{t) is also shown. 
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FIG. 2. The square amplitudes of the singlet states [Si) 
and 1 52) for the same situation as in figure |l|, but with a four 
times smaller ramp time of only r = Ti/Uh- The inset shows 
\ip2{t)f on a magnified scale. The dynamics of the system is 
clearly in the non-adiabatic regime. 
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FIG. 3. A square root of a swap, which is obtained from 
the situation of figure hi by halfing the pulse duration T. The 
probability of double occupancies is again strongly suppressed 
after the tunneling pulse. The resulting state is a fully entan- 
gled complex linear combination of \Si) and |r°), or, equiva- 
lently, of the incoming state ( |l9| ) and the outgoing state ( |2C| ) 
of the full swap. The quantum mechanical weigths of the 
latter states are plotted as thick solid lines 
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FIG. 4. A swapping processes in the absence of inversion 
symmetry along the axis connecting the dots. The square 
amplitudes of the singlet states \Si), i £ {1, 2, 3}, are denoted 
by I ^i I . The additional matrix elements entering the Hamil- 
tonian ^ ) are X = 0.2Uh and F = OAUh- The parameters 
of the tunneling pulse are r = Sh/Un and A = Uh /S with 
a duration T appropriate for a single swapping. As a result, 
a clean swapping operation can be performed also in the ab- 
sence of inversion symmetry. 



